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ID ■ 1 Introduction 

O 

Let O be a bounded region in and let H be a fourth-order, self-adjoint, uniformly 
elliptic operator on L 2 (Q) subject to Dirichlet boundary conditions on dQ, 



Hu{x) = D^{a Tl( (x)D <: u(x)} , 
i<?l<2 

ICI<2 



- T— I | 

The scope of this paper is to obtain integral boundary decay estimates for solutions of 
5-j \ the equation 

Hu = f, feL 2 (n). (1) 
More precisely, we want to establish ranges of f3 > for which the integrals 



j d~ 2 ~ fi u 2 dx and J d~ P \V 



u\ 2 dx 



(d(x) = dist(x, dfl)) are finite. If the boundary 0, is regular in the sense that the 
Hardy-Rellich inequality 

J^(Au) 2 dx>cJ^{^^ + ^}dx, ueH 2 {n), 

is valid, we then immediately have such an estimate since H$(Q) = Dom( J H" 1 / 2 ). Our 
aim is to establish better decay estimates that exploit the fact that the solution u of 
(prj belongs not only in Hq(O) but also in Dom(ff). 
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This problem is well studied in the case of second-order operators. In |D2| Davies 
obtained boundary decay estimates of the form 

+ ^p^dx <c(^\Hu\\ 2 \\H l / 2 u\\ 2 + \\u\\i^, ueDom(iJ), (2) 

for a > in some interval (0,ao). Here ceo is an explicitly given constant which de- 
pends on the boundary regularity and the ellipticity constants of H. As an application 
of ([2|) stability estimates were obtained on the eigenvalues {A n } of H under small 
perturbations of the boundary dQ. More precisely, if Q C £1 is a domain such that 
<9f2 C {x € £1 : dist(x, d£l) < e} and if {A n } are the corresponding Dirichlet eigen- 
values (the operator H being defined by form restriction), then it was shown that ([2]) 
implies 

< A n - A n < c n e 2a (3) 

for all n £ N and all e > small enough. This estimate has obvious applications in the 
numerical computation of eigenvalues; see |D2] for more on specific examples. 
Inequality © was subsequently improved in [D4], where e 2a , a < ao, was replaced by 
e 2a °, for the same Qo; this was done by estimating the integrals Jd(z)<e \ Vu\ 2 dx and 
Id(x)<e u2 dx for small e > 0. For results analogous to those of |D4] for the p-Laplacian 
together with applications we refer to Fleckinger et al. |FHT] . See also |EHK] where 
estimates of this type were first obtained for eigenfunctions of second-order operators. 
For relevant results in the case of singular operators see [M] . 

In our main theorem we obtain integral decay estimates analogous to ([2]) for fourth- 
order operators. More precisely, for a fourth-order operator H with L°° coefficients we 
establish boundary decay estimates of the form 

L {^ + ^ + ^) dx ^i^^ «€Dam(20, (4) 

for a in an interval (0, ao). Under additional assumptions we obtain ao = 1/2, which is 
optimal. To prove (pE]) we first use some general inequalities, which lead to a property 
(P a ) being identified as sufficient for the validity of Q. We then study property (P Q ), 
and find sufficient conditions under which it is valid; the distance function used here is 
taken to be the Finsler distance induced by the operator. 

Technical reasons oblige us to make a regularity assumption that is not needed in the 
second-order case and requires d(x) to be C 2 near d£l. This relates to a recurrent and 
largely unsolved issue in higher-order problems: a distance function is normally only 
once differentiable, but is required for technical reasons to be differentiated more than 
once; see for example [Bj , where such an issue has arisen in the context of heat kernel 
estimates. 

Finally, as an application of Q we obtain stability bounds analogous to ([3]) on the 
eigenvalues of H under small boundary perturbations. 

The structure of the paper is as follows: in Section [2] we provide a sufficient condition 
(P a ) for the validity of (j3J); in Section [3] we establish the range of a for which (P a ) is 
valid for different classes of operators; and in Section 2] we present the application to 
eigenvalue stability. 
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Setting 



We fix some notation. Given a multi-index rj = (rjx, . . . , r]]\f) we write r/\ = r/i! . . . rj^l 
and \rj\ = r)\ + • • • + t]n- We write 7 < rj if 7^ < rji for all i, in which case we also 
set d< = rj\/~f\{ri — 7)!. We use the standard notation D v u = (d/dxi)^ 1 . . . (d/dxN) r,N u 



C3 however are the same throughout the paper. 

We now describe our setting. We assume that £1 is a bounded domain in with 
boundary dQ. We consider a distance function d(-, •) on $7, and denote by d(-) the 
corresponding distance to the boundary d£l. We say that d(-) belongs in the class T> if 
it satisfies: 

(Dl) There exist c%, 02 > such that for any x,y £ £1 

(a) ci < \V z d(z,y)\ < C2 , z e $7 

(b) cid Eu c{x,y) < d(x,y) < c 2 d Eu c{x,y). 

(D2) There exist 6, r > such that 

(a) d(x) is C 2 on {0 < d{x) < 6} 

(b) \V 2 d(x)\ < cd(x)- 1+T on {0 < d{x) < 9}. 

(D3) The following Hardy-Rellich inequalities are valid for all v £ C 2 (£l): 



The distance cf(-, •) will typically be a Finsler distance, in which case (a) and (b) of 
(Dl) are equivalent. Condition (D2) is essentially a strong regularity assumption on 
dQ, as will be seen below. Its validity in examples will always involve the specific value 
r = 1; we choose however this more general and somewhat axiomatic setting because, 
we believe, it shows more clearly what the essential ingredients are. Finally, for more 
on Hardy-Rellich inequalities, optimal constants as well as improved versions of such 
inequalities we refer to |BFT[ IBT] and references therein. 

In the sequel we shall often need to twice differentiate d{x) near <9Q. In order to avoid 
repeatedly splitting integrals in two, we redefine d(x) on {x G 0, : d(x) > 6} so that 
now d(x) is C 2 and positive function on O such that d(x) equals inf{d(x, y) : y S dQ,} 
for x £ {dist(x, d£l) < 8} but not necessarily for all x G fl (of course d(x, •) extends 
to dO, by uniform continuity). In relation to this we emphasize that throughout the 
paper what really matters is what happens near the boundary d£l. We also note that, 
while the validity of estimate and assumption (D3) is independent of the specific 
distance d(-) £ T> chosen, we shall need to consider non-Euclidean distances since some 
of the intermediate calculations do depend on the specific choice of distance. 



and (Vit) ?? = u%\ . . .u™. By V 2 n we denote the vector (u XiXj )fj =1 . The letter c will 
denote a constant whose value may change from line to line; the constants ci,C2 and 




We note that a sufficient condition for (D3) (b) is the Hardy inequality 
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We will consider operators of the form 

Hu{x) = D ri {a v( (x)D' : u(x)} , x GO, (5) 

M=2 
ICI=2 

subject to Dirichlet boundary conditions on d£l; lower-order terms can be easily ac- 
comodated. More precisely, we start with a matrix-valued function a{x) = {a r) ^(x)} 
which is assumed to be have entries in L°°(f2) and to take its values in the set of all 
real, N(N + l)/2 x N(N + l)/2 matrices (N/(N + l)/2 is the number of multi-indices 
7] of length \rj\ = 2). We assume that {a r) ^(x)} is symmetric for all x G fi and define a 
quadratic form Q(-) on the Sobolev space Hq(Q) by 

Q(u)= [ Y a vC {x)D v u(x)D i u(x)dx, u G fln(fi). 

| T/|=2 

Kl=2 

We make the ellipticity assumption that there exist A, A > such that 

XQo(u) < Q(u) < AQ (u) , u € F 2 (O) , 

where Qo(u) = J n (Au) 2 dx denotes the quadratic form corresponding to the bilaplacian 
A 2 . We then define H to be the associated self-adjoint operator on L 2 (J7), so that 
(Hu,u) = Q(u) for all u G Dom(ff). 



2 Boundary decay 

Let d(-) G T>. Let a > be fixed and let us define 

u(x) = d(x)~ a , 

We regularize oj defining 

d n (x) = d(x) + - , = d n (x)~~ a , n = l,2, .... (6) 

n 

We note that u G Hq(Q) implies w ra u G Hq(£1), n G N. It is crucial for the estimates 
which follow that, while they contain the functions d n and u n , they involve constants 
that are independent of n G N. 

Lemma 1 Let a > 0. There exists a constant c which is independent of n G N suc/i 
that 

f I |V 2 u| 2 |Vd 2 

] fl h2^ + ^ + ^ H^ C ^M- ueH (n). (7) 



Proof. It suffices to prove (|7|) for all u G C 2 (f2). So let u G C 2 (0) be given and let 
v = uj n u, a function also in C 2 (0). Using (D3) we have 

-dx = I ^-rdx 



V 2 



n 



< c [ (Avfdx . 
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Similarly, 

I" iViil 2 f 1 

J n $+^ dx = J n -^\ud*- l vVd n + d*Vv\ 2 dx 

f v 2 , f \Vv\ 2 , 

< C / -rrdx + C / ^—dx 

Jn dl Jn d l n 

< c [ (Avfdx, 

Jn 

where we have used the fact that 

/ \V 2 v\ 2 dx= [ (Avfdx. (8) 
Jn Jn 

Finally, since d and d n differ by a constant, 

= d"v XiXj + ad°~ 1 (d Xi v Xj + d Xj v Xi ) + ad^~ x d XiXi v + a(a - l)d%~ 2 d Xi d Xj v, 
and therefore 

f \V 2 U\ 2 , f f ,„ 2 |2 J /■ |V«| 2 , /■ V 2 , /• |V 2 d| 2 2 , 1 

/ — i^—dx < c< \V v\ dx+ — ^— dx + / -rrcfe + / — ^ — 1> dx >. 
7n d 2a Un Jn dl Jn d± Jn d? n J 

Since d n > d, the second and third terms in the brackets are smaller than c f^{Av) 2 dx 
by the Hardy-Rellich inequalities (D3). The same is true for the last term by (D2). 
Thus, one more application of (JSj) concludes the proof. // 

Lemma 2 Let a G (0, 1) and uj n = d~ a . Then there exists a constant c > 0, indepen- 
dent of n G N, such that 

Q{u,oj 2 n u) < c\\Hu\\ 2 \\H a / 2 u\\2 , u G Dom(ff). 



Proof. For any n£N and u G C*c(^) we bave 

Jn 

= In^ dX 
< cQ(u). 

Hence tof/ a < c£f in the quadratic form sense, which by [Dl|, Lemma 4.20] implies that 
u>n < off* (since a G (0, 1)). Hence given u G Dom(ff) we have 

Q^tc^lt) < 11^11211^^112 

< c|| J ff«|| 2 ||-H' af/2 «||2, 

which is the stated inequality. // 
We can now establish a sufficient condition for the boundary decay estimates. 
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Theorem 3 Let a G (0, 1) be fixed and let uj n = d n a . Assume that there exist k, k' > 
independent of n G N such that 

Q{uj n u) < kQ{u,uj 2 n u) + k'\\u\\ 2 2 , u G C 2 (n), (9) 

for all n G N. T/ien i/iere exists c > sitc/i i/icrf 

^ + yfe + ^J^< c ll^ll2ll^ /2 ^ll 2) ^Dom(^). (10) 

Proof. The validity of Q for all u G C 2 (Q,) implies its validity for all u G Hq(Q) and 
in particular for all u £ Dom(ff). Hence given u G Dom(ff) and applying Lemmas Q] 
and [2] we conclude that there exists a constant c such that for any n G N there holds 

Letting n — > +oo, applying the Dominated Convergence Theorem and using the fact 
that the spectrum of H is bounded away from zero we obtain (fT0j) . / / 

3 The property (P a ) 

The validity of assumption ([9]) of Theorem [3] will be our main interest in this section. 
For the sake of simplicity, for any a G (0, 1) we define the property (P Q ) (relative to 
the distance function d G T>) as 

There exists constants k, k! > such that 

Q(d- a u) < kQ(u, d~ 2a u) + k'\\u\\l (P a ) 
for all n G N and u G Cl(Q). 

This is precisely assumption of Theorem [3l Our aim in this section is to obtain 
sufficient conditions under which property (P Q ) is valid. In the following three subsec- 
tions we present three theorems that provide such conditions. The first applies to all 
operators in the class under consideration; the second applies to operators of a specific 
type but gives a better range of a > 0; and the third applies to small perturbations of 
operators in the second class. 

Remark. If dQ. is smooth then the ground state <p °f A 2 decays as d{x) 2 as x — > d£l. 
Hence the integral in the left-hand side of (|10p is not finite for a > 1/2. For this reason 
and throughout the rest of the paper we restrict our attention to a G (0, 1/2). 

3.1 General operators 

We always work in the context described at the begining of Section [2j We recall that 
for a G (0, 1/2) we have uj n = d~ a = (d + ^)~ a ', we also recall that A and A are the 
ellipticity constants of the operator H. 
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Theorem 4 There exists a computable constant c > such that property (P a ) relative 
to the Euclidean distance is valid for H for all a G (0, c~ 1 A~ 1 \). 

Proof. Let u G C 2 (£l) be fixed. Setting v = d~ a u and using Leibniz' rule we have 

Q(d- a u) -Q(u,d~ 2a u) 

= Q{v)-Q{d a n v,d~ a v) 

= / E a vC (D^vD^v - D^{d^v)D^{d- a 

M=2 

ICI=2 

= " / E E c^4a v ^d a n )(D 5 d- a )(D^v)(D<- s v)dx 

M=2 1<V 
|CI =2 5<C 
7+<5>0 

< cA / E l yi(i nl • l vi O • |V 2 ~^| • \V 2 ~ j v\dx. 

0<i,j<2 
i+j>0 

But, by (Dl) and (D2), 



|Vd^ Q | = adf' 1 , \V 2 dt a \ < cadt a ~ 2 , 
and we thus obtain 

Q(v)-Q(dy,d- a v) < cAaJ^\V 2 v\ 2 + ^ + ^y x 

< cAX^aQiv). 

Hence, if a is such that cAA _1 a < 1, then property (P Q ) is valid for H. / / 

3.2 Regular coefficients 

The weak point of Theorem 0] is the poor information it provides on the range of a for 
which (P a ) is valid. In this subsection we shall consider operators of a more specific 
type and for which we shall see that (P a ) is valid for all a G (0, 1/2). 

It will be useful in this subsection to drop the multi-index notation and write the 
quadratic form as 

N 

Q(u) 



/ E a ijkiu XiXj u XkXl dx , u e Hq(Q) . 

n i,j,k,l=l 



We may clearly assume that the functions have the following symmetries: 

O'ijkl — Ojjikl > O'ijkl — ^ijlk j &ijkl — Qklij • (1-Q 

We make the following additional assumptions on the coefficients {aijki}' 

(i) There exist 9,t > such that 

(a) each a^ki is differentiable in {d(x) < 9} (12) 

(b) iVoyjwl < ccT 1+T on {d{x) < 9} 

(ii) E a ij ki{x)£,i£ k rij'qi < E ^jkli^^mVl > t V £ , x Eft . 
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Without any loss of generality we assume that r in (i) is the same as in (D2). Condition 
(ii) is a technical one, whose necessity is not clear. We present two examples in which 
it is valid. 

Example 1. Suppose that aijki = bijbki for some non- negative N x N matrix {bij}ij. 
Then (ii) is valid by the Cauchy-Schwarz inequality for the non-negative form (£, rf) i— > 
bij^ifjj. This for example includes operators of the form Aa(x)A, for which we have 
a-ijkl = a(x)5ij8ki- 

Example 2. Suppose that = ^i^M^/b where = > for i, k = 1,...,N. 
Then it is easily seen that (ii) is again valid. 

We choose the distance function d(-) to be the one naturally associated with H, that 
is the one induced by the Finsler metric p(x, rf) whose dual metric (cf. (|15p below) is 

i,j,k,l 

This implies in particular that the function d(-) satisfies 

Y a ijki(x)d Xi d Xj d Xk d Xl = 1 , a.e. iGil. (14) 

i,j,k,l 

Indeed, the inequality J2%,j,k,l a ijki{x)d Xi d Xj d Xk d Xl < 1 is shown in [A, Lemma 1.3]. To 
prove the reverse let y denote a point of differentiability of d. Then y has a unique 
nearest point yo G d£l; so d(y) = d(y,yo) =■ s. Let yt, t G [0, s], be the geodesic joining 
yo and y parametrised by arc length so that y s = y. Then for e > small we have on 
the one hand 

d(y s -e) ~ d(y) = d{y, y s - e ) = p(y, y - y s _ e ) + o(e), 
and on the other hand, by differentiability, 

d(y s - e ) ~ d(y) = Vd(y) ■ (y s - e - y) + o(e) . 

Hence 

P*(y, Vd y = sup — - > hm — = 1 . (15) 

p(y,0 p{y,y-ys-e) 

We note that the metric is Riemannian if the symbol of the operator H is the square 
of a polynomial of degree two. 

We assume that our basic hypotheses (Dl) — (D3) of the Introduction are valid; Con- 
cerning in particular the validity of condition (D2), we note that it is satisfied if enough 
regularity is imposed on the boundary and the coefficients. If for example the boundary 
is and the coefficients aijki lie in C 3 (0 < d(x) < 9), then d £ C 2 (0 < d(x) < 6); 
see [LN] Section 1.3]. On the other hand, for the Euclidean distance a C 2 boundary is 
enough [GT] P 354]. 

It is useful to introduce at this point a class A of integrals that are in a sense negligible. 

Definition. A family of quadratic integral forms T n (v), v G C 2 (J2), n G N, belongs to 
the class A if for any e > there exists c e > (independent of n G N) such that 

\T n (v)\ < eQ{v) + c € f v 2 dx , n G N, v G C 2 {n). (16) 
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Lemma 5 Let I n {v) = J b n (D" / v)(D s v)dx, \j\,\S\ < 2, be a term that results after 
expanding Q(d^v, d~ a v) and integrating by parts a number of times. If b n contains as 
a factor either a derivative of or a second-order derivative of d n , then (I n ) n £ A. 



Proof. After expanding Q(d%v, d~ a v ) (cf. (|19j) below) we obtain a linear combination 
of integrals, and direct observation shows that each one of them has one of the following 
three forms (we switch temporarily to multi- index notation): 



(a) J a vC d-^ +8 \ {Vd n y+^- 8 (D^v){D & v)dx , \ v \ = \(\ = 2 , 7 < r, , 5 < C : 

(b) / a vC d- 3 +W (Vd n )^(D^d n )v(D^v)dx , \r,\ = |C| = 2 , 7 < r, , 

(c) / a v(: d- 2 (D^d ri )(D^d n )v 2 dx , \rj\ = \(\ = 2. 



(These are distinguished by the number of second-order derivatives of d n that they 
contain - none, one and two respectively.) Hence all resulting integrals have the form 

b n (x)(D^v)(D s v)dx , 0<| 7 |,|5|<2, 

where b n (x) is a product of a„£ with powers and/or derivatives of d n and, since Vd n is 
bounded, 

K(x)\ <cd n (x)- i+ ^, If!!. (17) 

In cases (6) and (c) however, where b n (x) contains as a factor at least one second- 
order derivative of d n , it follows from condition (D2) of the Introduction that we have 
something more, namely 

K(x)\ < cd n (x)- 4+ l 7+<5 l +r , x G n . 

This easily implies that the integral lies in A in this case. 

Suppose now that we integrate by parts in the integral above, transfering one derivative 
from, say, D^v, (j 7 | > 1), to the remaining functions. If the derivative being transfered 
is d/dxi, we obtain - in an obvious notation - the integral 

{a r , c d- 4+ ^(Vd n ) T > + ^- s (D s v)} x (D^v)dx. 

If the derivative d/dxi "hits" either a„t or one of the factors that make up (V(i n ) ,?+ ^~ 7_5 
we obtain an integral of the form 

b n (x)(D^v)(D 5 v)dx 

JSl 

where \b n \ < cd~ 1+T (f n 4+ ' 7+<5 ' ; hence this integral belongs in A. // 

Example. We illustrate the last lemma with an example: in (I21j) below there appears 
the integral 

I n (v) = / a ijk id~ 2 d Xi d X:j vv XkXl dx . 
J o 
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Letting (T n ) n , (T^) n denote elements in A we compute 



a-ijkid n 2 d Xi d x , vv XkXl dx 

n 



(aijki)x k d n 2 d Xi d Xj vv Xl dx + 2 / aij k id n 3 d Xi d Xj d Xk vv Xl dx 



^•ijkid n d x . x .d x .vv x ,dx / &ijkid n d Xi d x . x ,vv x , xl dx 
n Jn 



n 



aijkid n 2 d Xi d Xj v Xk v Xl dx 



aijkid n 3 d Xi d X] d Xk (v 2 ) Xl dx - / ®ijkid n d Xi d x - v x v x dx + T n (v) 
n Jn 

= 3/ a ijM d~ A d Xi d X] d Xk d Xl v 2 dx - a ijk id~ 2 d Xi d Xj v Xk v Xl dx + T^(v) . 

Note. The summation convention over repeated indices will be used from now on. 

Lemma 6 There exists (T n ) n £ A such that 

Q(v) - Q(d%v,d~ a v) = 2a 2 / a ijM d~ 2 d Xt d Xi v Xk v Xl dx (18) 

Jn 

+4a 2 / a ijk id~ 2 d Xi d Xk v x .v Xl dx - (a 4 + 11a 2 ) / d~ A v 2 dx + T n (v) , 
Jn Jn 

for alive C 2 (Cl). 

Proof. For (3 equal to a or to -a we have 

{d P n v) XiXj = d^+pdt'd^+pdt'd^+piP-^d^d^v + pdt'd^v . 

(19) 

We substitute in Q(d"v, d n a v) and expand. Now, by Lemma [5] all terms containing 
second-order derivatives of d n belong to A. Further, the symmetries (fTT|) of give 

Q j ijkid X id X jd Xk v Xl — Q'ijkid Xi d Xk d Xl v X j — ... 

Q'ijkid Xi d X jV XkXl — Q,ij k id Xk d Xl v XiX j — ... (20) 

Q-ijkid Xi d Xl v X jV Xk — Q>ij k id Xi d Xk v X jV Xl — .... 

Denoting by (T n ) n an element of A which may change within the proof we thus arrive 
at 

Q(d*v,d~ a v) = J^a ijk i^v XiXj v XkXl + 2a 2 d~ 2 d Xi d X] vv XkXl - Wd~ 2 d Xi d Xk v Xj v Xl 

+<ia 2 d~ 3 d Xi d X] d Xk vv Xl + (21) 
a 2 (a 2 - l)d~ A d Xi d Xj d Xk d Xl v 2 Xdx + T n (v). 

We integrate by parts the second and fourth terms in the last integral. By Lemma 
[5j all terms that contain either derivatives of a{j k [ or second-order derivatives of d n , 
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belong in A. Hence, denoting always by (T n ) n a generic element of A we obtain (cf. 
the example above) 

aijkidn 2 d Xi d Xj vv XkXl dx = 3 / a ijk id~ 4: d Xi d Xj d Xk d Xl v 2 dx 

- / a ijk id~ 2 d Xi d x .v Xk v Xl dx + T n (v) 
Jn 

and, similarly, 

/ a ijk id~ 3 d Xt d Xj d Xk vv Xl dx = - j a ijkl d~ 4 d Xt d Xj d Xk d Xl v 2 dx + T n (v). 

Jn ^ Jn 

Substituting in (j2~Tj) yields 

Q(d%v,d~ a v) = J^a ijk i^v XiX] v XkXl - Wd~ 2 d Xi d Xk v Xj v Xl 

-2a 2 d~ 2 d Xi d Xj v Xk v Xl + (a 4 + lla 2 )d~ 4 d Xi d Xj d Xk d Xl v 2 jdx + T n (v). 
Recalling that ((111), relation {18]) follows. // 

Lemma 7 Lei u E C 2 (f2) and u; = d n 3 ^ 2 v. Then there holds 



Q(v) - Q(d%v, d n a v) = 2a a ijkl d n d Xi d Xj w Xk w Xl dx 

+4a 2 [ a ijM d n d Xi d x k w Xi w Xl dx + (-a 4 + -^-) / d~ 1 w 2 dx + T„,(u), 
where (T n ) n £ -4. 

Remark 1. When working with the function u>, all integrals have the form 

/ b n (x)(D^w)(D 5 w)dx 
Jn 

where the function b n satisfies 

K{x)\ < cd n {x)- l+ ^ +s \ (22) 

Such an integral lies in A if in addition there holds 

Mar) | < cd n (x)" 1+l7+<5|+T , x E . 

for some r > 0; as before, these are precisely the integrals that contain either second- 
order derivatives of d n or (first-order) derivatives of a,ij k i. 

Remark 2. Since d and d n differ by a constant, for the sake of simplicity we shall write 
d Xi instead of (d n ) Xi , etc. 

3 i 

Proof of Lemma^ We substitute v Xi = d%w Xi + ^dnd Xi w in Recalling the symme- 
try relations (|20p (with w in the place of v) and using the fact that aij k id Xi d Xj d Xk d Xl = 1 
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we obtain 



Q(v)-Q(d«v,d- a v) 
4a 2 I a ijld d~ 2 d Xi d Xk 



n 



3 31 1 r - 3 - 1 

dn w Xj + -dn d Xj w dnW Xl + -dnd Xl w dx + T n (v) 



+2a 2 / a ijk id n 2 d Xi d Xi 
Jn 



dnW Xk + -d%d Xk w 



dnW Xl + -dhd Xl w 



dx 



-(a 4 + lla 2 ) / d~ 1 w 2 dx + T n (v) 
Jn 



4a l^aijkid n d x% d Xk d n w Xj w Xl + 3d n d Xj w Xl w + -d n d Xj d Xl w z 
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dx 



+2a / a ijM d n d Xi d x 



dlw Xk w Xl + 3d 2 n d Xk w Xl w + -d n d x ,d x ,w 2 



dx 



-(a 4 + lla 2 ) / d~ l w 2 dx + T n (v) 
Jn 

= 4a 2 / aij kl d n d Xi d Xk w Xj w Xl dx + 2a 2 / aijkid n d Xi d Xj w Xk w Xl dx 

+ (-" 4 + -^-) f d~ 1 w 2 dx + 18a 2 / a ijk id Xi d Xi d Xk w Xl wdx + T n (v). 
2 Jn Jn 

But the last integral belongs in A by an integration by parts; hence the proof is com- 
plete. // 

Lemma 8 Let v £ C 2 {Q) and w = d n 3 ^ 2 v. Then there holds 



Q(v) 



aijkid\w XiX] w XkXl dx + - / aijkid n d Xi d Xj w Xk w Xl dx 



2Jn 
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-3 / a ijk id n d Xi d Xk w Xj w Xl dx + / d n 1 w 2 dx + T n (u) , 
/n Id 



where {T n ) n is an element of A. 



Proof. We have 



v 



3 3 ± 3i 3-i 3i 

d*w XiX . + -d%d Xi w Xj + -dld X] w x% + -d n 2 d Xi d X] w + -d%d XiXj w 



As already mentioned, all terms involving second-order derivatives of d n belong in .4.. 
Hence, using the symmetry relations (|20p once more we compute, 



Q(v) = I aijki 



'3 1 3-i 

dnW XiXj + 3dnd Xi w Xj + -d n 2 d Xl d Xj w 

3 1 3-1 

dnW XkXl + 3dnd Xk w Xl + 2 d Xk d X[ w 



dx + T n (u ) 



(23) 



Oyfci d\w XiX] w XkXl + 6d 2 d Xi w Xj w XkXl + ^d n d Xt d Xj w XkXl w 

9 9 

+9d n d Xi d Xk w Xj w Xl + -d Xi d X] d Xk w Xl w + — d~ d Xx d Xj d Xk d Xl w z 



dx + T n (i>) 
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The fifth term belongs in A be a simple integration by parts. We also integrate by 
parts the second and third terms, obtaining respectively 



aijkid n d Xi w Xj w XkXl dx = -2 \^a ijk id n d Xi d Xk w Xj w Xl dx 

+ a ijk id n d Xi d Xj w Xk w Xl dx + T n (v) 



aijkld n d Xl d Xj w XkXl wdx = - I a ijk id n d Xi d Xj w Xk w Xl dx + T n (v) . 

Substituting in (|23j) and recalling (fl~4"|) we obtain the stated relation. / / 

We can now prove the main theorem of this subsection. For any a £ (0, 1/2) we define 



k 9 



l-4a 2 )(9-4a 2 )' 



Theorem 9 For the operator H and relative to the metric i f 13\) property (P a ) is valid 
for all a £ (0,1/2). More precisely, for any a £ (0,1/2) and any k > k a there exists 
k! < +oo such that 

Q(d~ a u) <kQ(u,d~ 2a u) + k f \\u\\ 2 2 , u£C 2 (n). (24) 



Proof. Let u € C 2 (f2) be given and let v and w be defined by v = d~ a u and w = d n 3 ^ 2 v 
respectively. Define j a = (40a 2 — 16a 4 )/9 and observe that j a £ (0)1)- Applying 
Lemmas [7] and [U and assumption (|12h (ii)we obtain 

la Q(d- a u) - [Q(d- a u) - Q(u,d- 2a u)\ 

= j a Q(v)-[Q(v)-QKv,d- a v)] 

= la J a ijk idnW XtXj w XkXl dx + - 2a 2 ) J a ijk id n d Xi d X:j w Xk w X[ dx 

-(37 Q + 4a 2 ) J a ijkl d n d Xi d Xk w Xj w Xl dx (25) 

+ (a 4 - + -j^j a ijk id~ x d Xi d X] d Xk d Xl w 1 dx + T n (v) 
> 7a / a ijk id z n w XtXi w XkXl dx + (^y- - 6a 2 ) ( a ijk id n d Xi d Xj w Xk w Xl dx 
+ (° 4 ~~ + "Iff) X a ijkidn 1 d Xt d Xj d Xk d Xl w 2 dx + T n (v) . 

Therefore 

7a Q(d" a n) - [Q(d~ a u) - Q(u, d~ 2a u)} > T n {v) , (26) 

since the coefficient of the last integral is zero and those of the other two integrals are 
non-negative. Now, for any e±, e 2 > we have from (|16p 

\T n (v)\ < eiQ(v) + c ei \\v\\l 

= eiQ(^n" n ) + c ei\\dn a u\\l 

< eiQ(d~ a u) + c ei (e 2 \\d~ a ~ 2 u\\ 2 2 + c e2 \\u\\ 2 2 ) 

< eiQ(d~ a u) + c ei (ce 2 Q{d~ a u) + c £2 \\u\\fj , 
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and therefore 

\T n (v)\ < eQ(d- a u) + c £ \\u\\ 2 2 (27) 
for any e > small. Choosing e > so that 7« + e < 1 we obtain from (|26[) and (I27p 

Q(d- a u) < - Q(u, d- 2a u) + - ^ . 

1 - 7 a - e 1 - 7 a — e 

Hence (|24p is valid with A; any number larger than 1/(1 — 7 Q ) = k a . / / 

3.3 Small perturbations 

In this subsection we prove a stability theorem on the validity of (P„). We denote by 
M + the cone of all coefficients matrices for the operators under consideration, that is 



M + = |a = {a^} \ v \ = \^ = 2 ■ Offc symmetric, real valued and measurable 

with XQ (u) < Q(u) < \()u(ii). a € C'r(O) (A. A > 0; 

equipped with the uniform norm 

Halloo := esssup|a(x)|oo ; (28) 

here |a(x)| is the norm of the matrix a(x) = {a v ((x)} v ^ considered as an operator 
on R iV ( Ar + 1 )/ 2 . We recall that A, A, etc, denote the lower ellipticity constant for the 
operators induced by the matrices a, a, etc. We have 

Lemma 10 There exists a computable constant c > such that for all a € (0,1/2) 
there holds 

" |V 2 (d»| • \V 2 {d~ a v)\dx < cQ (v) , v G C 2 (n). 



Note. For an estimate on the constant c see the remark at the end of this subsection. 
Proof. For any (3 E R we have 

{d?v) 9iXj = d^v^+pdt'd^+pdt'd^+PiP-^dt'd^v + pdt'd^v . 

(29) 

We write this for (3 = a and for (3 = — a, and we multiply the two relations; cancels 
with d~ a and we obtain 

,^2/ 1 m ,^2/ m f,^2 ,2 IVvl 2 V 2 |V 2 d n | 2 o] 

|v 2 K»l • |v 2 K^)| < c I |v 2 H 2 + i-^L + ^ + L -^r L « 2 1 • 

The proof is concluded by using assumption (D2) on V 2 d and the Hardy-Rellich in- 
equalities (D3); here we have also used the fact that J Q \V 2 v\ 2 dx = f n (Av) 2 dx. // 
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Proposition 11 Let a G (0, 1/2) be fixed. Assume that (P tt ) is valid for the matrix 
a G M + relative to some distance d(-) G 2? and Ze£ k,k' > be such that 

Q(uj n u) < kQ(u,oj^u) + k'\\u\\l , reGN, ii£C c 2 (0). 

T/ien i/iere exists a constant c > swc/i i/tai i/d G M + satisfies \\a— a||oo < A[(l+cA;)] _1 , 
i/ien (P Q ) is also satisfied for a relative to d(-); more precisely, there exists k' < +oo 
so that 

k 

Q(cj n u)<- =— - — —Q{u,u 2 n u) + k'\\u\\l, neN, u G C 2 (f2). 

1 - A -i (l + cfc)||a - a||oo 



Proof. We first note that 

\Q(v) -Q(v)\ < [ \a- a\-\V 2 v\ 2 dx <\\d- aW^Qoiv) , v£C 2 (n). (30) 

Moreover, setting v = u) n u we have from Lemma [TUI 

IQCm,^ 2 ^) — Q(u,u 2 u)\ < \\a — alloc / |V 2 uj • |V 2 (w 2 -u)|dx 

= ||d- alloc / (V 2 ^ 1 ?;)! • \V 2 (uj n v)\dx 
Jn 

< c\\a - a\\ooQo(v) . (31) 
Prom (I30j) and (|3ip we conclude that for any n G N and u G C 2 (fi) we have 

Q(uj n u) < Q(uj u u) + \\a - a\\ooQo(uj n u) 

< kQ(u,uj 2 u) + fe'||M||| + ||d - a||oo<3o(^n^) 

< k(^Q(u,u; 2 u) + c|[a - aWocQoiiOnu)^ + A:'||n|| 2 + ||a - a\\ooQo{u n u) 

= kQ(u,uj 2 u) + (1 + ck)\\a - a||oo<3o(w n n) + fc'IMI 2 

< kQ(u,uj 2 u) + A~ 1 (l + c/c)||d — a\\ooQ(uj n u) + /c'||tt||2, 

from which the statement of the lemma follows. / / 

Let Q denote the cone of all coefficient matrices that satisfy assumptions (i), (ii) of 
Section [3.2i Also let k a be as in Theorem [9j Combining Proposition [11] and Theorem 
Owe obtain immediately the following 

Theorem 12 There exists a computable constant c > such that if for some a G 
(0, 1/2) the coefficient matrix a of the operator H satisfies 

dist^oo(a, Q) < 



(1 + ck c 
then (P a ) is satisfied for H . 

Proof. Let a G Q be such that 

A 
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By Theorem [9] (P Q ) is satisfied for a and (|24p is valid for any k > k a . If in addition k 
satisfies 

A 

l|a_a||o ° < TTck' 

then (P Q ) is also valid for a by Proposition (TTJ / / 

Example. Suppose that the coefficients a„£ are uniformly continuous and satisfy 
(|12p (ii). Then property (P a ) is valid for H for all a £ (0,1/2). This is seen by 
approximating a„£ with smooth functions using an approximate identity; note that the 
approximating functions also satisfy (ii). 

Remark. The constant c of the above proposition is precisely the constant c of Lemma 
[TU1 Precise estimates for this constant can be easily obtained. Indeed, it follows from 
(1291) that for a G (0, 1/2) there holds modulo A 

72/JQi„.M IT72/ 



V\d%v)\-\V\d-<*v)\dx 

Vv l2 .,2 



< 3/ (\V 2 v\ 2 + 4 a 2 \Vd\ 2 ^ + a 2 (a + l) 2 \Vd\^ 

Jn v a* d\ 

< 3 /(|W + |v,|^4lv,|^) & 



2 9 „ v 2 

( iv^r + |vc<r 

m v 

Hence, letting 02,03 be as in (Dl), (D3), we obtain (modulo .4) 

jT |V 2 (d»| • |V 2 (d-%)|dx < 3(l + cleg 1 + ^cg 1 ) J^(Av) 2 dx . 

In fact, since we work modulo 4, the last constant can be improved 3(1 + c\A~ x + 
(9/16)c2-E> -1 ), where A and B are the weak Hardy constants, that is they satisfy 

iVvpdx > A I -prdx — c v 2 dx 
Jn d A Jn 



[ (Avfdx > B [ V -rdx - c" [ v 2 dx . 
Jn Jn « 4 Jn 



For smooth boundaries with a smooth Riemannian metric this amounts to A = 1/4, 
B = 9/16. 



4 An application: eigenvalue stability 

In this final section we demonstrate how the boundary decay estimate of Theorem 
[3] yield stability bounds on the eigenvalues {A n } of H under small perturbations of 
the boundary d£l. The proof follows closely the corresponding proof in [D2| for the 
second-order case, however we include it here for the sake of completeness. So we 
consider a distance function d(-) £ T>, an operator H as above and assume that the 
boundary decay estimates (|10p are valid for some fixed a € (0,1/2). For e > we 
define Q e = {x G £1 : d(x) > e}. We assume that e < 6/2 so that d(x) is C 2 on Q\Q,2e- 
We define d e (x) = dist(x, c?0 £ ), x € 0, and make the additional assumption that there 
exists c > such that for small enough e > there holds 

\V 2 d t \ < c on {d(x) < 2e}. (32) 
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Now, let Q be a domain such that 

Q e C C ; 

we do not make any regularity assumptions on dQ. We denote by {A n } the eigenvalues 
of the operator H on L 2 (f2), which is defined by restricting the quadratic form Q(-) on 

h 2 (q). 

Let <?!> be a non-negative, smooth, increasing function on R such that 



m = | 

We define a C 2 cut-off function r on 11 by 

t(x) = 



0, t < 0, 

1, t>l. 



o, x £n\n e , 

4>(d e (x)/e), x £ Q e . 



Note that t(x) = 1 when > 2e; moreover (j32[) yields 

|r(x)|<l , |Vt(x)| < ce" 1 , |V 2 r(x)| < ce~ 2 . 

Let us now denote by {4> n } the normalized eigenfunctions of H. For n > 1 we set 

L n = lin{0i, . . .,0 n } , L„ = lin{r0i, . . . ,T0 n }, 
and observe that L n C H 2 ((l). 

Lemma 13 There exists a constant c > such that for all e > sma// and all u G 
Dom(if) i/iere ZioWs 

(i) |Q(ru) - Q(u)| < ce 2a \\Hu\\ 2 \\H a l 2 u\\ 2 , 



(ii) 



T!l b - It 2 



Kce^WHuW^WH^uW 1 ^ 



Proof. Let u € Dom(i/). On f^e we have t(x) = 1, hence 



\Q(tu)-Q(u)\ 



2 {p"(r«))(D f (ru)) - (Z>"u)(Z> c u)} da; 

|7?|=2 
ICI=2 

< c / ||V 2 (ru)| 2 + |V 2 u| 2 )<£c 

Jd(x)<2e 1 J 

< c/ i|V 2 -u| 2 + |Vr| 2 |Vn| 2 + |V 2 r[V)cfe 

Jd(x)<2e L > 

< c[ {\y 2 u\ 2 + e~ 2 |Vn| 2 + e~ 4 u 2 \dx 

Jd{x)<2e 1 > 

|V 2 d 2 



< ce 



.{ 



|Vu| \u\ 



x)<2 6 >- d 2a d 2+2a d i+2a 



17 



from which (i) follows by means of Theorem [3l Similarly, 



2 <r II 112 

< \\TU — U\\2 



< / |w| 2 (ix 

Jd(x)<2e 

A+2a I u<1 



Jd(x)<2e d 4+2a 

< ce 4+2a \\Hu\\ 2 \\H a / 2 u\\ 2 , 
from which (ii) follows. // 

Theorem 14 Assume that there exists a distance function d G T> and an a G (0, 1/2) 
such that M0\) is satisfied. Assume also that \313il is valid. Then there exists c, d > 
such that for each n > 1 there holds 

< A n < A n < A n + cA^ 4 e 2Q , (33) 
for all e > satisfying e 2a < c'Xn 5 ^. 

Proof. We fix n > 1. Since L n C -Hq(^) we nave by min-max 

A n < sup{Q(tO/|M|2 : v G L n } 

= sup{(5(™)/||rn||2 : u G L n }. (34) 
Now, let u G L n be given. It follows from Lemma [13] (i) that 

Q(tu) < Q(u) + ce 2a \\Hu\\ 2 \\H a / 2 u\\2 

< Q(n) + Ce 2Q A n / 4 ||n|| 2 . (35) 

Similarly Lemma PT3l (ii) gives 

||rn|| 2 > || U || 2 - C e 2+a A n / 8 || U ||2(||«||2 + ||™||2) 

> || U || 2 - Ce 2+ -A n / 8 || U || 2 . (36) 

Assumimg in addition that ||u||2 = 1 we thus obtain from (|35p and (|36|) that 

Q(tu) Q{u) + ce 2 »A n /4 

I ~ 1 _2+ a \V8 



I"*ll2 l-ce 2 +°A 



< Q(«) + cA n / 4 e 2a 

< \ n + c\^e 2a , 

where for the second inequality we have used the fact that e 2a < c'X n 5 ^, with d small 
enough but fixed (and independent of n and e). Hence (|34p implies 



A„ < A n < A n + c\ 5 J* 



:2a 



which completes the proof of the theorem. / / 

Remark. In the case where Cl = B(l) and £1 = B{\ — e) we have X n = (1 — e) _4 A n 
and hence A n — A n = 4A n (e + 0(e 2 )). Hence the value a = 1/2 is the best possible for 
estimate (1331). 
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